Two-scale competition in phase separation with shear 
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The behavior of a phase separating binary mixture in uni- 
form shear flow is investigated by numerical simulations and 
in a renormalization group (RG) approach. Results show 
the simultaneous existence of domains of two characteristic 
scales. Stretching and cooperative ruptures of the network 
produce a rich interplay where the recurrent prevalence of 
thick and thin domains determines log-time periodic oscilla- 
tions. A power law growth R{t) ^ t" of the average domain 
size, with a = 4/3 and a = 1/3 in the flow and shear direction 
respectively, is shown to be obeyed. 

PACS numbers: 47.20Hw; 05.70Ln; 83.50Ax 



The application of a shear flow to a disordered binary 
mixture quenched into a coexistence region greatly af- 
fects the phase-separation process . A large anisotropy 
is observed in typical patterns of domains which appear 
greatly elongated in the direction of the flow |^ . This be- 
havior has consequences for the rheological properties of 
the mixture. A rapid strain-induced thickening followed 
by a gradual thinning regime is observed Various 
numerical simulations confirm these observations ^-|^. 

In a recent paper 0| we have studied the phase- 
separation kinetics in the context of a self-consistent ap- 
proximation, also known as large-iV limit. Within this 
approach the existence of a scaling regime characterized 
by an anisotropic power-law growth of the average size of 
domains was established. The segregation process, how- 
ever, cannot be fully described by this technique because 
interfaces are absent for large- A^ 

The aim of this letter is to investigate the ordering 
process by extensive numerical simulations. Our main 
result concerns the simultaneous existence of two length 
scales which characterize the thickness of the growing 
domains. The competition between these scales produces 
a rich dynamical pattern with an oscillatory behavior due 
to the cyclical prevalence of one of the two lengths. A 
power law growth R{t) ~ of the average domain size, 
with a — 4/3 and a — 1/3 in the flow and shear direction 
respectively, is shown to be obeyed by a renormalization 
group (RG) analysis. 

The kinetic behaviour of the binary mixture is de- 
scribed by the Langevin equation 



(1) 



where the scalar field tp represents the concentration dif- 
ference between the two components of the mixture . 
The equilibrium free-energy can be chosen as usual to be 
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where b, k > and a < in the ordered phase, v is an 
external velocity field describing plane shear flow with 
average profile given by 



(3) 



where 7 is the shear rate and Cx is the unit vector in 
the flow direction. 77 is a gaussian white noise, represent- 
ing thermal fluctuations, with mean zero and correlation 
{T]{r,t)ri{r',t')) = -2TTV^5{r- r')5(t - t'), where F is 
a mobility coefficient, T is the temperature of the fluid, 
and the symbol (...) denotes the ensemble average. 

The Langevin equation (|^) has been simulated in d — 
2 by first-order Euler discretization scheme. Periodic 
boundary conditions have been used in the flow direc- 
tion while in the y direction the point at {x, y) is iden- 
tified with the point at {x + jL,y + L), where L is the 
size of the lattice [§. Lattices with L = 1024, 2048, 4096 
and space discretization intervals Ax — 0.5, 1 were used. 
In the phase separation process the initial configura- 
tion of iy9 is a high temperature disordered state and 
the evolution of the system is studied in model (0) with 
a < 0. Parameterization invariance of (|l|) allows one to 
set F=|a|=6 = K=l. The structure factor is de- 
fined as C{k,t) = {Lp{k,t)Lp{—k,t)) where <p{k,t) are the 
Fourier components of ip. Results will be shown for the 
case 7 = 0.0488, T = 0, L = 4096, Ax = 1, (ip) = 0. 
Similar results have been obtained for other values of the 
parameters. 

A sequence of configurations at different values of the 
strain jt is shown in Fig. 1. After an early time, when 
well defined domains are forming, the usual bicontinu- 
ous structure of phase separating domains starts to be 
distorted for jt^l. 
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jt = 20 





FIG. 1. Configurations of a portion of 512 x 512 sizes of the 
whole lattice are shown at different values of the strain 'yt. 

The growth is faster in the flow direction and domains 
assume the typical striplike shape aligned at an angle 0(t) 
with the direction of the shear which decreases with time. 
As the elongation of the domains increases, nonuniformi- 
ties appear in the system: Regions with domains of dif- 
ferent thickness can be clearly observed at jt = 11. The 
evolution at still larger values of the strain is shown for 
'yt = 20. The domains with the smallest thickness even- 
tually break up and burst with the formation of small 
bubbles. 

A systematic existence of two scales in the size distri- 
bution of domains is suggested by the behavior of the 
structure factor, shown in Fig. 2. At the beginning (see 
the picture at — 0.2) C{k, t) exhibits an almost circu- 
lar shape, corresponding to the early-time regime without 
sharp interfaces. Then shear-induced anisotropy becomes 
evident, C(fc, t) is deformed into an ellipse, changing also 
its profile and, for 7/: > 1, four peaks can be clearly ob- 
served. The position of each peak identifies a couple of 
typical lengths, one in the flow and the other in the shear 
direction. The peaks are related by the k ^ —k symme- 
try so that, for each direction, there are two physical 
lengths. This corresponds to the observation of domains 
with two characteristic thicknesses, made in Fig. 1. 

The dips in the profile of C(/c, t) develop with time 
until C{k,t) results to be separated in two distinct foils 
at 'jt ~ 4. The evolution of the system imtil this stage 



is well described by the solution of the linear part of 
(letting 6 = 0): 



C(fc, t) — Co exp" 



- k^(s){k^(s)-a)ds 



(4) 



where k(s) — {kx,ky + ^sk^) and Co is the structure 
factor at the initial time Then non- linear effects 

become essential in producing the patterns shown in 
Fig. 2 at 7t — 11,20. At = 11 we evaluate the 
positions of the peaks at {kx,ky) — (0.015,0.107) and 
{kx,ky) ~ (0.038,0.23). This gives a value around two 
for the ratio between the characteristic sizes of domains; 
the same value is found for the ratio between the posi- 
tions of the two peaks in the hystogram of the domain 
size distribution. The relative height of the peaks in one 
of the foils of C{k,t) can be more clearly seen in Fig. 3 
where the two maxima are observed to dominate alter- 
natively at the times jt = 11 and = 20. 
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FIG. 2. The structure factor is shown at different values of 
the strain 'yt. 

The competition between two kinds of domains is a 
cooperative phenomenon. In a situation like that at 

— 11, the peak with the larger ky dominates, describ- 
ing a prevalence of stretched thin domains. When the 
strain becomes larger, a cascade of ruptures occurs in 
those regions of the network where the stress is higher 
and elastic energy is released. At this point the thick 
domains, which have not yet been broken, prevale and 
the other peak of C{k,t) dominates, as at jt = 20. In 
the large- limit the prevalence of one or the other peak 
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has been shown to continue periodically in time . Here 
we have an indication of a similar behavior, although the 
observation of the recurrent dominance of the peaks on 
longer timescales is hardly accessible numerically. 



7i = 11 



7^ = 20 



FIG. 3. The same structure factor of Fig. 2 at "ft 
in a three-dimensional plot. 
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The hallmarks of this dynamics are found in the 
behavior of the average size of domains, Rx(t) and 
Ry{t), in the flow and shear direction. In our sim- 
ulations these quantities have been calculated through 
R^{t) = {J dkklC{k,t)/ J dkCik^t))-"^/^, and analo- 
gously for Ry{t); their behavior is shown in Fig. 4. Due 
to the alternative dominance of the peaks of C{k,t), 
and Ry increase with an oscillatory pattern. The lat- 
ter reaches a local maximum at the characteristic times 
tm when C{k,t) is of the form of Fig. 2 at 7t = 20 
and thick domains are more abundant. The time in- 
terval T„ — tn+i — tn between two cascades of rup- 
tures increases exponentially as the phase separation pro- 
ceeds so that the oscillations appear to be periodic on 
a logarithmic time scale. This is expected because do- 
mains are growing and longer and longer times are re- 
quired to break them. The following argument sim- 
ply illustrates the origin of the log-time periodic oscil- 
lations: The break up of an elongated domain is caused 
by the shear flow, which enters Eq. (|l|) through the last 
term on the l.h.s., whose magnitude we infcre to be pro- 
portional to T~^. For a sharp interface exposed with 
an angle 0{t) to the flow this term is proportional to 
7sin0(i) ~ lO{t) ^ ^ Ry{t) / Rx{t) , where the asymptotic 
smallness of 9{t) has been taken into account. The RG ar- 
gument developed below shows that ^Ry{t)/Rx{t) ~ 
so that T„ ~ tn and hence lnt„ ^ n + cost.: The rupture 
events occur periodically in h\t. 

Stretching of domains requires work against surface 
tension and results in an increase A?] of the viscosity 
[^,|o|. We calculate the excess viscosity as Ari{t) — 
-7"^/ j0^kxkyC{k,t) 0. Starting from zero At? 
grows oscillating up to global maximum at 7^ ~ 12. Then 
the excess viscosity relaxes to zero in a similar way to 
what observed in previous simulations Q]. The relative 
maxima of A77 are found in correspondence of the minima 
of Ry when the domains are maximally stretched. 

In the case without shear the asymptotic kinetics is 



characterized by dynamical scale invariance, which is re- 
flected by a power law growth R[t) ~ of the average 
size R[t) of domains ||]. The value of the growth ex- 
ponent a is related to the physical mechanism operating 
in the separation process and is a = 1/3 for diffusive 
growth. In the case with shear the existence of a similar 
behavior has not been clearly assessed. A stationary state 
characterized by domains of finite thickness is generally 
observed in experiments [ pj| . Only in some experimen- 
tal realizations the existence of a regime with power law 
growth has been shown . 

The self-consistent solution of Eq.(|l|) shows Q that a 
generalized scaling symmetry holds when a shear flow is 
applied and different exponents are found for the power- 
law growth of Rx{t) and Ry{t). However these exponents 
are correct for models with vectorial order parameter and 
do not directly apply to the case of a binary mixture Q . 
The existence of a scaling symmetry and the actual value 
of the growth exponents cannot be reliably studied by 
direct numerical simulation of Eq.(|l|): This is because, 
even for the bigger lattice we have used, the fast growth 
in the z-direction makes finite size effects relevant be- 
fore a full realization of the scaling regime occurs. More- 
over, the oscillatory behavior of Rx{t) and Ry{t) prevents 
a straightforward computation of the growth exponents, 
unless the evolution of the system is followed over sev- 
eral decades, which is, presently, impossible. Therefore, 
in order to infer the actual value of the growth exponents 
we resort to a RG analysis. 

As usual, we define the RG transformation for the 
Fourier components '^-^(t) — 1/y/V J df(p{f,t)e'^'''^ of the 
field 95 (r, t). They verify the equation 
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with < f]k{t)rik'{t') >= 2Trk^S{t ~ t')d{k + k'). Taking 
into account the anisotropic growth of domains, we gen- 
eralize the RG scheme of by considering the change 
of scale in 



kx ^ k^ — kxb ^ , 

t^t' 

and the field transformation 



ky — kyb° 



tb- 



^~,{t) = b^v'~(t') 



(6) 



(7) 



where b is the rescaling factor and the meaning of the 
exponents ax^OLy^C, will be clarified in the following. Di- 
mensional analysis implies that the structure factor can 
be written as 



C{k,t)^Rx{t)Ry{t)f{x,y) 



(8) 



with X = kxRx(t),y = kyRy(t). Its invariance in form 
with respect to the transformations (0,^) gives C = {oix + 
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ay) 12. The invariance of the scahng variables x^y under 
the transformations (|^) imphes that Ux , a.y correspond to 
the growth exponents. 
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FIG. 4. Evolution of the average domain size in the shear 
(lower curve) and flow (upper curve) directions. The straight 
line has slope 4/3. 



a RG approach. Results show the simultaneous existence 
of domains of two characteristic sizes in each direction. 
The two kind of domains alternatively prevail, because 
the thicker are thinned by the strain and the thinner are 
thickened after cascades of ruptures in the network. This 
mechanism produces an oscillation which decorates the 
expected power-law growth R{t) ~ of the average size 
of the domains, with a = 4/3 and 1/3 in the flow and 
shear direction, respectively. The oscillations occur on 
logarithmic time-scale as in models describing propaga- 
tion of fractures in materials where the releasing of elastic 
energy is measured Finally, in a recent paper [ p^ , 
a first treatment of the effects of hydrodynamics on the 
phase separation of a binary mixture in uniform shear 
has been given. It would be interesting to know if and 
how hydrodynamics affects the global picture described 
in this letter. 
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The dynamical scaling regime is described by a fixed 
point of the Langevin equation under the set of trans- 
formations (^,0) ijl^. Generalizing the prescription of 
[ p^ , we assume that the fixed point hamiltonian T veri- 
fies the relation 



(9) 



meaning that the main contribution to the free energy is 
due to the interfaces in the direction of the flow. By 
inserting (^,0) into (||) and using (||), one obtains a 
Langevin equation similar in form to (^ with rescaled 
parameters 



T' T 



T' = Tb- 



(10) 
(11) 



A fixed point of the above recursions with F, 7 7^ is 
obtained when ay — 1/3, Ux — ay + 1, with the temper- 
ature being not relevant for the process of phase separa- 
tion. We observe that a difference between the growth 
exponents Aa ~ ax — ay in the range 0.8 1 has been 
measured in ||l^,|l^. Moreover, the above analysis sug- 
gests that, for a constant value of the strain jt, the excess 
viscosity scales as A?/ ~ 7"^ with /? = 1/3. 

In conclusion, we have studied the phase-separation ki- 
netics of a binary fluid in an uniform shear flow by direct 
numerical simulation of the constitutive equations and in 
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